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Abstract. In this paper wc prove that the center of mass of a supercritical 
branching-Brownian motion, or that of a supercritical super-Brownian motion 
tends to a limiting position almost surely, which, in a sense complements a 
result of Tribe on the final behavior of a critical super-Brownian motion. This 
is shown to be true also for a model where branching Brownian motion is 
modified by attraction/repulsion between particles. 

We then put this observation together with the description of the interact- 
ing system as viewed from its center of mass, and get the following asymptotic 
behavior: the system asymptotically becomes a branching Ornstein Uhlenbeck 
process (inward for attraction and outward for repulsion), but 

(1) the origin is shifted to a random point which has normal distribution, 
and 

(2) the Ornstein Uhlenbeck particles are not independent but constitute a 
system with a degree of freedom which is less by their number by precisely 
one. 



1. Introduction 

We start with some basic notation. 

Notation 1. In this paper MfiR'^) andA^i(R'^) denote the space of finite measures 
and the space of probabihty measures, respectively, on WK For n e MfiW^), we 
define ||/i|| := pi{M.'^). 

1.1. A model with self-interaction. Consider a dyadic (i.e. precisely two off- 
spring replaces the parent) branching Brownian motion (BBM) in M.'^ with unit 
time branching and with the following interaction between particles: if Z denotes 
the process and Zl is the ith particle, then Zl feels the drift 

- E A^' - 

l<.j<nt 

where 7^0, that is the particle's infinitesimal generator is 
(1.1) + i E 7(^.^-)-V. 

l<J<"t 

(Here and in the sequel, nt is a shorthand for 2 L*J , where [t\ is the integer part of 
t.) If 7 > 0, then this means attraction, if 7 < 0, then it means repulsion. 
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To be a bit more precise, we can define the process by induction as follows. Zq 
is a single particle at the origin. In the time interval [m, m + I) we define a system 
of 2™ interacting diffusions, by the following system of SDE's: 

Zl = W^ + ^ r (^:-^:)ds, re [0,1); * = 1,2,...,2- 

l<j<2" 

where {W^ ,t £ [0, 1)} are independent Brownian motions for i = 1, 2, . . . , 2™ and 
they start at the position of their parents at the end of the previous step (at time 
m — 0). Existence and uniqueness follows from the fact that these 2™ interacting 
diffusions can be considered as a single 2™rf-dimensional Brownian motion with 
drift b : M^'"'' K^^d ^^^^j.^ 

b(a;i,X2, ...,Xd,xi+d,X2+d, ■■■,X2d, a;i+(2'"-i)d, 2:2+(2"-i)d, •••,a;2™rf) 

= : (/3i,/32,...,/32"d)^ 

satisfies 

(3^ = 2"™(xi + Xi+d + •■• + 2:i+(2'"-i)d) - Xi, 

(32 = 2~"\X2 + X2+d + ■■■ + X2+{2"--l)d) - X2, 
Pd = ^^"^{Xd + X2d + ■■■ + X2'^d) - Xd, 

f3i+d = 2"™(a;i + xi+d + ■■■ + a;i+(2"«-i)<i) 

/32+d = 2""(.X2 + X2+d + ■■■ + X2+{2'"-l)d) 
P2d 2^"^{xd + X2d + ■■■ + 2^2™ d) - X2d, 

f3i+{2"--i)d = 2""(a;i + xi+d + ■■■ + 2^i+(2™ 

P2+{2'"-l)d = 2""(a;2 + X2+d + ■■■ + X2+(2"- 
/32'^d = 2"™(2:<i + X2d + ••■ + X2'^d) - X2^d- 

Existence and uniqueness is guaranteed by the Lipschitzness (linearity) of b. 

Remark 2. It may seem natural to replace the interaction we defined by the 
gravitational force between particle^, however this would lead to a randomized 
version of the (notoriously difficult) 'n-body problem.' O 

1.2. Results on the self-interacting model. We are interested in the long time 
behavior of Z , and also, whether we can say something about the number of particles 
in a given compact set for n large. In the sequel we will use the standard notation 

{Zug) ■■^T.Z.gizi). 

In this paper wc will show (Theorem [T4)) that Z asymptotically becomes a 
branching Ornstein Uhlenbeck process (inward for attraction and outward for re- 
pulsion), but 

(1) the origin is shifted to a random point which has d-dimensional normal 
distribution 7V''(0, 2), and 



- Xi+d, 

- X2+d, 



-i)d) ^ a;i_|_(2m_i)d, 

-l)d) - X2+{2'"-l)d, 



I.e. the forces are given by Newton's law of universal gravitation: they vary as the inverse 
square of the distance between the particles. 
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(2) the Ornstein Uhlenbeck particles are not independent but constitute a sys- 
tem with a degree of freedom which is less by their number by precisely 
one. 

For the local behavior we formulate and motivate a conjecture (Conjecture I15|). 

1.3. An extension of Tribe's result on critical super-Brownian motion. In 

the proof of Theorem [M] we will first show that Zt '■— -^J27=i^t' the center of 
mass for Z satisfies hmt^oo Zt = N, where N ~ JV^{0,2). In fact, the proof will 
reveal that Z moves like a Brownian motion, which is nevertheless slowed down 
tending to a final limiting location (see Lemma [5] and its proof). 

Since this is also true for 7 = (BBM with unit time branching and no self- 
interaction), our first natural question is whether we can prove a similar result for 
the supercritical super-Brownian motion. 

Another motivation for the same goal is as follows. Tribe [7J proved that a critical 
super-Brownian motion near its extinction time ^ behaves like a single Brownian 
path stopped at ^. More precisely, Aj ^ <5f as i — > 00 a.s. in the weak topology, 
where is a d-dimensional random variable and its distribution is the same as that 
of a Brownian motion at time ^. 

We would like to extend Tribe's result to the supercritical super-Brownian motion 
A. We are interested in whether we can obtain a similar result on the survival set. 
Of course, then Xf does not shrink to a point in any sense, however, we may hope 
to get an analogous result regarding the center of mass, defined as 



(Since f{x) = x is not a bounded function we may not hope to use Tribe's techniques 
though.) 

In the sequel we will prove that A will be a time changed Brownian motion on 
a finite (but random) time interval. 

At the beginning of this subsection we referred to a result on Z. In fact we now 
see that the results on A and Z are completely analogous: in both cases the center 
of mass is a Brownian motion slowed down in such a way that the time interval 
[0, 00) is compressed into a finite one. A slight difference is that, in case of Z, 
the terminal time is deterministic, because so is the offspring distribution. (The 
terminal time is t — 2.) 

Let A be the (iA, /?, a; M'') -superdiffusion with a, /? > (supercritical super- 
Brownian motion). Let P denote the corresponding probability. Let us restrict H. 
to the survival set 



Since (3 > 0, F^(5) > for aU fi^O. 

Our main result is that the center of mass for A stabilizes as t —^ 00, and 
furthermore the path of the center of mass is a finite piece of a Brownian path 
(with a different time parametrization) . 

Theorem 3. Let a, /? > and let X denote the center of mass process for the 
P, a; M."^)- superdiffusion X . Then, onS, 

(i) Af converges Ps^-almost surely as t ^ 00. 

(ii) In fact, the finite path {Af}t>o is the same as the path of a time changed 
Brownian motion on [0,Too), where Too is a positive and finite random 




5 G 17 I Af(cj) > 0, > 0}. 
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variable. More precisely, there exists a d- dimensional Brownian motion B 
(on a,n enlarged space) such that X = B o T where T : M_|_ — > is a 
random time change satisfying limt^oo T{t) < oo Ps^ — a.s. 

Remark 4. 

(a) A heuristic argument for (i) is as follows. Obviously, the center of mass 
is invariant under ff-transforms whenever H is spatially (but not tem- 
porarily) constant. Let H{t) := e-f^K Then is a (i A, 0, e^'^^'a; K'')- 
superdiffusion, that is, a critical super-Brownian motion with a clock that 
is slowing down. Therefore, heuristically it seems plausible that , the 
center of mass for the transformed process stabilizes, because this is obvi- 
ously true in case of extinction, and otherwise the center of mass under the 
heat flow does not move. 

(b) From (ii) one can easily conclude some of the a.s. path properties of X. 
For example, since the p-variation (j> > 0) is invariant under changing the 
parametrization, we know that any segment of the path has infinite total 
variation, but the whole path has finite quadratic variation a.s. <> 

2. The mass center stabilizes 

Notice that 



- E N-^*) = -f E 



(2.1) - >^ [zi-zi) = -\ zl-mzi\=z,-zi 

^l<j<nt 

and so the net attraction pulls the particle towards the center of mass (net repulsion 
pushes it away from the center of mass). Thus the following lemma is relevant: 

Lemma 5 (Mass center stabilizes). Let Zt ■= 5^"=i ^^^.^ ^ is the center 
of mass for Z. Then limt^oo Zt = N, where N ~ A/''^(0, 1). 

Proof. For t € [m, m + 1), there are 2™ particles moving around. U t,t + At € 
[m,m+ 1), then Zl moves as a Brownian motion plus a vector 

^(Zt - Zi)At, 

so 

Zl+At = Zl + i(Zt - Zl)At + Bi,{At), 
where the {Bl{s), s > 0; i = 1,2,. ..,2™} axe independent Brownian motions 
starting at the origin. Thus 

2" / 2™ \ 2™ 

Zt+At = ^t+At = Zt+ 72-" 2"Zt -J2zi]At + 2-" i3^(At), 

i=l V i=l J i=l 

where denotes independent sum. Since 2"'Zt - Ei=i = 0, we obtain that for 
< T < 1, 

2m 

(2.2) Zm+r = Zm + 2"™ 0B^(r). 

Let B'^"'\t) := 2-"/^ 0^^"^ B^(t). Using induction, we obtain that 

% = mm^m\l)®. . . ^BC'Hi) ■ ■ •e^=|=B(^^"-^)(l)®-^B(-)(l), 
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(it is easy to check that, as the notation suggests, the summands are independent) 
which, using Kolmogorov's Three Series Theorem converges almost surely. (We 
will denote the a.s. limit by N.) On the other hand, since _B(™) is a Brownian 
motion, we can apply Brownian scaling and get 2~'"/^i?('"^ (r) = W^"^^ ( 2^) ' where 
are independent Brownian motions. We have 



and so, 



b(")(i) e ^B(i)(i) ® • • • ® — i — ® ^b("')(i) 

V2 V2L*J-i ^/ii7 



where t := t — nt. Since the summands are independent on the right hand side, N 
has the same distribution as a Brownian motion a,t t — 2, that is, N ^ Af{0, 2). □ 

For another proof see the remark after Lemma [S] 

Remark 6. It is interesting to note that Z is in fact a Markov process. Indeed, 
the distribution of Zt {t = m + t, < t < 1) conditional on J>„ is the same as 
conditional on Zm, because Z itself is a Markov process. But the distribution of 
Zt only depends on Z,„ through Zm, as 

Z.iz„®iy(^'")(^), 

whatever Zm is. O 

We will also need the following fact later. 
Lemma 7. The coordinate processes of Z are independent. 

We leave the simple proof to the reader. 

3. Normality via decomposition 

We will need the following result. The decomposition appearing in the proof will 
also be useful. Recall that m = [ij . 

Lemma 8. {Z^, Zf,...,Zf ) is joint normal for all t > 0. 

Proof of lemma. By Lemma [T51 we may assume that d = 1. We prove the 
statement by induction. 
For m = 1 it is trivial. 

Suppose that the statement is true for m — 1. S ince (^Zm—i^Zm—±^...TZm—i^ is 
normal, we can consider it just as well as a 2™ dimensional degenerate normal at 
the instant of the fission of the particles. Indeed, the vector 

/ r^l r^l ryl ry^ ry'2^^ \ 

l^r?j.-lJ ^m-lJ ^m-lJ ^m-li ^rn-l) 

has the same distribution on the 2™^^ dimensional subspace 

S := {x eM."^ \ xi = X2,X3 = X4, ...,X2'r^^i = 2:2™} 

of M^'" as the vector ^/2{Zm-i, Z^-i: Z^_i) on R^'" \ (The reader can eas- 
ily visualize this for m — 2: the distribution of [Z\,Z\) is clearly V2 times the 
distribution of a Brownian particle at time 1, i.e. AA(0, ^/2) on the line xi — X2.) 
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Since the convolution of normals is normal, therefore, by the Markov branching 
property, it is enough to prove the statement when the 2™ particles start at the 
origin and the clock is reset: t € [0, 1). 

Define the 2™ dimensional process Z* on the time interval t G [0,1) by 

starting at the origin. Because of the interaction between the particles attracts the 
particles towards the center of mass, Z* is a Brownian motion with drift 

j[(ZuZu...,Zt)-(ZlZl...,Zf 

Notice that this drift is orthogonal to the vectoJl v :— (1, 1, 1), that is, the vector 
{Zt, Zt, Zt) is nothing but the orthogonal projection of {Z}, Zf, Z^ ) to the 
line of V. This observation immediately leads to the following decomposition. The 
process Z* can be decomposed into two components: 

• the component in the direction of v is a Brownian motion 

• in the ortho-complement of v, it is an independent Ornstein-Uhlenbeck 
process with parameter 7. 

It follows from this decomposition that Z* is Gaussian. □ 

Remark 9. Consider the Brownian component in the decomposition appearing in 
the proof. Since, on the other hand, this coordinate is 2™/^Zf, using Brownian 
scaling, one obtains another way of seeing that Zt stabilizes at a position which is 
distributed as the time l-|-2~^-|-2~^-|-... + 2^™ + ... = 2 value of a Brownian motion. 
(The decomposition shows this for d = I and then it is immediately upgraded to 
general d by independence.) o 

Corollary 10 (Asymptotics for finite subsystem). Let k > 1 and consider the 
subsystem {Z^ , Zf, Z^), t > mo for mo := [logfcj + 1. (This means that at time 
mo we pick k particles and at every fission replace the parent particle by randomly 
picking one of its two descendants.) Let the real numbers ci,...,Cfc satisfy 

k k 

(3.1) ^c. = 0, ^c,2 = l. 

i=l 4=1 

Define xj^j'^i' - ''^'') ^^^^ '^i^t '^'^'^ note that is invariant under the translations 
of the coordinate system. Let Ct denote its law. 

For every k > 1 and ci, Ck satisfying 13. - is the same d- dimensional 
Ornstein Uhlenbeck process corresponding to the operator 1/2A — 7V-x, and in par- 
ticular, 

lim = A/- ( 0, — 

t^oo V 27 

For example, taking ci — l/\/2, C2 = ^ l/\/2, we obtain that when viewed from 
a tagged particle's position, any given other particle moves as \/2 times the above 
Ornstein Uhlenbeck process. 

Proof. By independence (Lemma I13p it is enough to consider d = 1. For m 
fixed, consider the decomposition appearing in the proof of Lemma|8]and recall the 



^For simplicity, we use row vectors in this proof. 
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notation. By p.l[) . whatever m > mo is, the 2™ dimensional unit vector 

(ci,C2, ...,Cfc,0,0, ...,0) 

is orthogonal to the 2™ dimensional vector v. This means that - is a one 
dimensional projection of the Ornstein Uhlenbeck component of Z*, and thus it is 
itself a one dimensional Ornstein Uhlenbeck process (with parameter 7) on the unit 
time interval. 

Now, although as m grows, the Ornstein Uhlenbeck components of Z* are defined 
on larger and larger spaces (S C K.^ is a 2™~^ dimensional linear subspace), 
the projection onto the direction of (ci, C2, c^, 0, 0, 0) is always the same one 
dimensional Ornstein Uhlenbeck process, i.e. the different unit time 'pieces' of 
V]y(ci,...,cfc) obtained by those projections may be concatenated. □ 



4. The interacting system as viewed from the center of mass 



Recall that by (|2.2p the interaction has no effect on the motion of Z. Let us see 
now how the interacting system looks like when viewed from Z. 

4.1. The description of a single particle. Using our usual notation, assume 
that t e [m,m + 1) and let r := t — [<J . Also recall that © denotes independent 
sum. When viewed from Z, the time r relocatioiH of a particle can be described as 
follows: 

2" 

A{Zl - Z,) = AZ} - AZ, = BI{t) - 2— BI{t) - ^{Z] - Z,)t. 

i=l 

So if Z^ - Z, then 

2"* 

AF/ = Blir) - 2-™ i?^(r) - jY,\. 

i=l 

Clearly, 

2" 2" 

Blir) - 2-"0B^(r) = 02-"B^(r) (1 - 2-")Bi(r), 

1=1 4 = 2 

and thus the right hand side is a Brownian motion with mean zero and variance 
d{l - 2-™)t. That is, denoting cr„ := 1 - 2"™, 

AYt^ ^darnW\T)-^Yt'T, 

where is a standard Brownian motion. 

We have thus obtained that (for fixed m) Y^ corresponds to the operatoiQ 

i(T„iA - 7x • V, 

which is an Ornstein-Uhlenbeck process on the time interval [m,m + 1). Since 
for TO large am is close to one, the relocation viewed from the center of mass is 
asymptotically governed by an 0-U process corresponding to ^A — 7a: • V. 



■^I.e. the relocation between time m and time t. 

^Here A is the Laplace operator unlike in previous lines where it denoted difFerence. 
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Remark 11 (Asymptotically vanishing correlation between driving BM's). Let 

= {W'^^ fc = 1, 2, d). For 1 < i 7^ j < 2", we have 



E 

2- 



B'o\r) - 2-™0BS(r) (r) - 2-"0BS(r) 



- 2" 



r=l 



■ 2"V = (2^-™ - 2-™)t = 2-"V, 



Var(Bf(r))+Var (i?f(r)) 
that is, for i ^ 

(4.1) E [W^^\T)W=-\r)\ = Sm ■ _ ^_„,)2 ^- 

Hence the pairwise correlation decays to zero as t ^ oo (recall that m = \t\ and 
r = t-me [0,1)). 

And of course, for the variances we have 

(4.2) E [W''''{t)W'^'^{t)] = Ske ■ (1 - 2"") r, for 1 < i < 2". o 

4.2. The description of the system; the 'degree of freedom'. Fix m > 1 and 

for t e [to, to + 1) let Yt :— {Yf^, Y^ )-^, where ()-^ denotes transposed. (This is 
a vector of length 2™ where each component itself is a cZ dimensional vector; one 
can actually look at it as a 2™ x d matrix too.) We then have 



where 



W\t) :=i3^(r)-2-™0B^(r), z = 1, 2, 2" 

are mean zero Brownian motions with correlation structure given by (j4.ip - (j4.2p . 

Just like at the end of subsection 11.11 we can consider F as a single 2™d- 
dimensional diffusion. Each of its components is an Ornstein-Uhlenbeck process 
with asymptotically unit diffusion coefhcient. 

By independence, it is enough to consider the d = \ case, and so from now on, 
in this subsection we assume that d —\. 

Let us first describe the distribution oiWt for t > fixed. Recall that {i?Q(s), s > 
0; i — 1,2,. ..,2™} are independent Brownian motions starting at the origin. By 
definition, Wt is a 2'"-dimensional multivariate normal: 



Wt 



/ 1-2- 



V 



_2-™ 
1 - 2-"' 



- 2"™ 



1 



Bo{t) A(")i?o(t) 



However, since we are viewing the system from the center of mass, it is a singular 
multivariate normal. Its 'true' dimension is the rank of the matrix A^"^\ 



Lemma 12. ran/t(A(")) = 2" - 1. 
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Proof. We will simply write A instead of A*^™). Since the columns of A add up 
to zero, the matrix A is not of full rank: r{A) < 2™ — 1. On the other hand, 

/I 1 ... 1 \ 



2™A 



1 1 



= 2™I, 



V 1 1 ... 1 / 

where I is the 2™-dimensional unit matrix, and so by subadditivity. 



(A) + 1 = r(2"A) + 1 > 2'^ 



□ 



By Lemma[T21 Wt is concentrated on the (2™ — l)-dimensional linear subspace given 
by the orthogonal complement of the vector (1, 1, 1)-^; in this 2™ — 1 dimensional 
subspace Wt has non-singular multivariate normal distribution. What this means 
is that even though Wi, W2, W^2™ are not independent, their 'degree of freedom' 
is 2™ — 1, i. c. the 2™-dimensional vector W is determined by 2™ — 1 independent 
components (corresponding to 2™ — 1 principal axes). 



5. Asymptotic behavior 

How can we put together that Zt tends to a random final position a.s. with the 
description of the system 'as viewed from Zj?' 

Lemma 13. For t >0, the random vector Yt is independent of the path {Zs}s>t- 

Proof. First, for any t > 0,Yt is independent of Zt, because (assuming d = 1) the 
vector 

i^t, z] — Zt, z1 — Zt, ■ ■ ■ ,z1 — Zt)^ 

is normal (since it is a linear transformation of the vector [Z\, Z^, . . . , Z^ Y' , which 
is normal by Lemma [5]), and so it is sufficient to show that Zt and Z'l — Zt are 
uncorrelated for 1 < z < 2™. But this is obvious, because the random variables 
Zl,Z'^ , Z'l are exchangeable and thus, denoting n = 2™, 

E\Zt{Z\-^t)\^E 

= U {{Z\f\ + ^E{ZlZl) [[Zlf] - !^^I^E{ZlZf) = 0. 

To complete the proof of the lemma, recall Lemma [5] and its proof and notice that 
the distribution of {Zs}s>t only depends on its starting point Zt, as it is that of a 
Brownian path appropriately slowed down, whatever Yt (or even Zt) is. Since, as 
we have seen, Yt is independent of Zt, we are done. □ 

We have the following result on the asymptotic behavior of the system. 

Theorem 14. The asymptotic behavior of Z is that of a branching Ornstein- 
Uhlenbeck process, but with the origin shifted by a random, normally distributed 
vector. More precisely, for almost all x, conditionally on limt_>oo Zt = x, one has 
Zt — Zt®Yt, where Y is the branching Ornstein-Uhlenbeck process of section lJj] 



Z^ 



zi + zt + ... + z- 
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Proof. We already know from Lemma[5]that the almost sure limit N := limt^oo Zt 
exists. Let 

P^(.) :=P(. \ N^x). 

By Lemma [T51 Yt is independent of N, and thus P^{Yt G •) = P{Yt E ■) for almost 
all X. By definition, Zf = Zt -\-Yt, and by the above discussion, we can in fact write 
Zt ~ Zt ® Yt under . Putting this together with the description of Yt in section 
14. li we are done. □ 
We close this section with a conjecture on the local behavior of the system. 

Conjecture 15. Let g e C+(R''). 

(i) // 7 > (attraction), then there exists a random variable N ^ Af'^{0,2), 
such that, conditional on N = xq. 



hm 2 



-"(Z„,.g) = ^(J)'^%xp(-7| --xop) 



Also, there exists a random variable M ^ H'^ ^0, ^2 + 1^^ (with cor- 
responding law V), such that 

\im 2-"E{Zn,g) ^E(M,g) a.s. 

n — ►CO 

(Here Id denotes the d-dimensional unit matrix.) 
(ii) 7/7 < (repulsion), then 

lim 2-"(Z„,5) = a.s. 

n — >oo 

Remark 16. The function 2 + in Conjecture fTST i') is monotone decreasing. 
The intuitive meaning is that stronger attraction results in smaller variance of the 
limiting distribution. o 

Explanation of the conjecture: Once we have Theorem UM we can try to put 

it together with the Strong Law of Large Numbers for the local mass from [1] for 
the process Y. So N is the final position of the center of mass. 

First, let 7 > 0. If the components of 1" were independent and the branching rate 
were exponential. Theorem 6 of [1] would be readily applicable. However, since the 
2"* components of Y are not independent (as we have seen, their degree of freedom 
is 2™ — 1) and since, unlike in [T], we now have unit time branching, the method of 
[1] must be adapted to our setting. This adaption would require some extra work. 

Once this is done, however, it immediately follows that there exists a random 
variable N ^ J\f'^{0, 2), such that, conditional on — xq, 

lim 2'"(Z„,5) = ^(J)''^'exp(-7|.-xon,<?^ a.s 



Thus the distribution of hm„^oo 2~"(Z„, g) is the convolution of A/'(0, 2/^) and 
AA(0, (27)"^/^), yielding the second statement in the first part. 

The 7 > case is similar but since the limiting density is translation invariant, 
i.e. Lebesgue (see Example 11 in [T]), the final position of the center of mass plays 
no role. 
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6. Proof of Theorem [3] 

(i) Since a, (3 are constant, the branching is independent of the motion, and 
therefore N defined by 

Nt := e-P'WXtW 

is a nonnegative martingale (positive on S) tending to a hmit almost surely. It 
is straightforward to check that it is uniformly bounded in and is therefore 
uniformly integrable (UI). Write 

' e-nXt\\ Nt ■ 

We now claim that A^oo > a.s. on S. Let A := {iVoo — 0}. Clearly SZ C A, and so 
if we show that P{A) = P(S'C), then we are done. As is well known, P(S'C) — e 
On the other hand, a standard martingale argument (see the argument after formula 
(20) in 2J) shows that < u{x) := — \ogPg^{A) must solve the equation 

^Aw + (3u — au^ ~ 0, 

but since Ps^{A) — P{A) constant, therefore — \ogPs^{A) solves (iu ~ av? = 0. 
Since N is UI, no mass is lost in the limit, giving P{A) < 1. So u > 0, which in 
turn implies that — logP^^ [A] — (3 /a. 

Once we know that N^o > a.s. on S , it is enough to focus on the term 
e~^*{x,Xt). Let H{t) e"^*. Then is a (iA, 0, e-^*a; E'*)-superdiffusion, 
that is, a critical super-Brownian motion with a clock that is slowing down. One 
can write 

e-''\x,Xt)^{x,X,"). 

Define Tg = := e^^^Ss; then the semigroup {Ts}s>o corresponds to Brownian 
motion. In particular then 

(6.1) T4id]=id, 

where id(a;) = x. Therefore (x, Xj^) is a martingaleQ If we show that the martingale 
is UI, we are done. It is enough to show that it is uniformly bounded in L^. To 
achieve this, define by gnix) — \x\ ■ l{|j,|<„}. Then we have 

EsA^,Xl')'^Es,\{x,X^)\' < EsMX,")\ 
and by the monotone convergence theorem we can continue with 

= lim £;5„(g„(x),Xf)2. 

n — *oo 

Since gn is compactly supported, there is no problem to use the moment formula 
and continue with 

ft \ ft 



lim (l+ [ dse-^'{So,T,[agl])) = 1 + a lim / ds e-P'T,[glm. 



Recall that {Ts;s > 0} is the Brownian semigroup, that is, Ts[f]{x) — Ex{Ws), 
where W is Brownian motion. Since gn{x) < \x\, therefore we can trivially upper 
estimate the last expression by 



1 + alim / dse->^'Eo{W^) = l + a [ ds se" 



^It does not matter that the function is unbounded and changes sign. 



12 



JANOS ENGLANDER, 



Since this upper estimate is independent of t, wc are done: 

t>o a 

(ii) Keeping in mind that we are on the survival set S := {u! G ft \ Xt{uj) > 

0, yt > 0}, we first claim that X is an (S*, P, {jFt}f>o)-martingale, where Tt ■— 
a{{Xt; t>0}). _ 

Let us start with showing that -El^tl < oo. To see this, first recall that 

ix,Xt)/\\Xt\\^^-^^. 

Thus, 

E[\Xt\;S] ^E[N^'\{x,Xi')\;S] , 

and since we have seen in part (i) that {x,Xp) is (uniformly) bounded in L^, by 
Cauchy-Schwartz, it is enough to see that 

E [iV-2; 5-] = ir^; 5-] < oo, Vt > 0, 

1. e. that S] < oo, Vt > 0. This is true because is a one-dimensional 
diffusion on [0, oo) with generator 2; (a ^^-1-/3^), which, on S, tends to infinity, and 
therefore, by Fatou's Lemma, limt_^oo S*] = 0. Hence 5] < 00 
for large fs, but then by continuity, S] < 00, Vi > 0. 

Next, since < 00, in fact EXt = for i > 0, because Xt is symmetrically 

distributed {Xt = —Xt). Indeed, 



\Xt\\ J \\\Xt\ 
where B* := {/x* \ ^ e B} and ^*(-) := 

We now show that E{Xt \ J^s) = for < s < i. By the Markov branching 
property, 

E(Xt I Ts) = Ej^^ J^^ Xts Xsidx), 

where X^'' are independent copies of super-Brownian motions starting at Sx ■ Since 
X is mean zero, shifting by x and using Fubini's Theorem, we can continue the last 
displayed formula with 

= llhi[ ^Xs{dx) = Xs. 

\\^s\\ 

Next, we show that X has continuous paths. 
Let 

u{x,t) := a;l{|^|<(^+g)t} + (^2^+ ^)t'^{\x\>{^+e}t}i 
where e > is fixed for the rest of the argument. By [5] it follows that there exists 
an a.s. finite random time T = T(lu) such that for all t > T, 

{x,Xt)^{u[-,t),Xt). 
(This is obviously true on the extinction set.) Hence, for t > T, 

{x, Xt) - {x, Xt+At) = {ui-,t), Xt) ~ (u(-, t + At), Xt+At) = 

{{u{-,t), Xt) - {u{; t), Xt+At)) + {{u{; t), Xt+At) - {u{; t + At), Xt+At)) 

/ + //. 
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Now / — + as At J, a.s. by the continuity of the paths (the continuity is in 
the weak topology and u{-,t) is a bounded continuous function for aU t > 0), and 
// — > as At I a.s., because 

//= (7(Ai),Xt+At), 
where |7(A<)| < {^/2l3 + e)At a.s. So 

// < {^ + e)At-\\Xt+At\\ a.s., 

and we are done since 

hin |lXf+Ai|| = ll^tll a.s., 

again because of the path continuity. 

This proves right continuity of the paths of X; the left continuity is similar. 

Once we know that the center of mass is a continuous martingale, we can utilize 
the multidimensional Dambis-Dubins-Schwarz Representation Theorem, which we 
cite below. (For the theorem and its proof see e.g. Theorem 18.4 in ^.) 

A d-dimensional continuous local J^-niartingale M = {M-^, AP, M'^) is said to 
be isotropic if [M^] — [M^] = ... = [M''] holds almost surely, where [M*] denotes 
the total variation process for AP. 

Proposition 17 (Theorem DDS). Let M be an isotropic continuous local T- 
martingale in M.'^ with Mq, and define 

T, := M{t > 0; [M^]t > s}, Gs J^t,, s> 0. 

Then there exists Brownian motion B (with respect to a standard extension 

ofG), such that a.s. B = MoT on [0, (M^)oo) and M = Bo [M^]. 

In order to finish the proof of Theorem [3](ii), first let d = 1. Since isotropy 
automatically holds, if M := {x,X") and := mi{t > 0;[M]t > s}, then by 
Theorem DDS, there exists a one-dimensional Brownian motion B (on an enlarged 
space) such that Z = B o T. 

Now let d >2. Note, that in our case 

{x,X,") = ((a:i,Xf),(x2,Xf),...,(xd,Xf)), 

if a; = {xi,X2, ...,Xd), and thus M :— {x,X^) is isotropic, because {xi,X^) 
is a time-changed Brownian motion (1 < i < d) by the already proven one- 
dimensional case and so [{xi, X-^)]t = Tt for all i a.s. Hence, if := {xi,X^) and 
Ts := mf{t > 0; [M^]t > s}, then by Theorem DDS, there exists a d-dimensional 
Brownian motion B (on an enlarged space) such that Z = B oT. 

Since we already know from part (i) that Zt has a finite a.s. limit, it follows that 

Too ■= lim Tt < oo, a.s. □ 

t — >oo 
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